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A NOTE ON A BROKEN LAYER 
IN AN ORTHOTROPIC LAMINATE COMPOSITE 

K. ARIN 

Lehigh University, Bethlehem, Pa., 18015, U.S.A. 


ABSTRACT 

An orthotropic laminate composite containing a completely 
broken layer is considered. The problem is formulated in 
terms of integral transforms and then reduced to a singular 
integral equation which is solved numerically. The strength 
of stress singularity at the crack tip is determined from a 
characteristic equation which is obtained by studying the dom- 
inant part of the singular integral equation near the end points. 
The stress intensity factors are given for various material 
properties . 

1. INTRODUCTION 

In a previous work by Arin [1] , the problem of an ortho- 
tropic laminate composite containing a layer with a crack and 
bonded to two half-planes of dissimilar materials has been 
solved and the effect of the material properties as well as 
various crack sizes on the stress intensity factor at the tip 
of the crack illustrated. It has been shown that there are 
basically two types of • orthotropic materials. The results were 
given for two different cases (Materials Type I and II) . It 
was concluded that the decrease in the stress intensity factors 
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was apparent in the case of a matrix of stiffer material. 

This is meaningful from the point of view of crack arrest. 

In this paper, the limiting case for which the crack 
reaches the interface will be investigated. It can be shown 
that the Fredholm kernels in [1] become unbounded as the crack 
touches the interface, hence requiring a separate treatment 
of the problem. Extracting the additional singularities, the 
singular integral equation of [1] can be solved in a similar 
manner. Here, only the material Type I will be considered. 

The solution for material Type II can be obtained similarly. 
The notation and certain results of [1] will be used wherever 
necessary. Both the plane strain and the generalized plane 
stress conditions are studied simultaneously. 

The isotropic counterpart of this problem has been solved 
by Gupta [2] and Ashbaugh [3] . 

2. FORMULATION OF THE PROBLEM AND THE SOLUTION 


The singular integral equation of the problem of an or- 
thotropic laminate composite containing a layer with a crack 
of length 2a is given as [1] : 


i r a ft ( t) dt + r a 
ir -a t-x -a 


k(x,t) ft (t)dt 


(1 " v xv v vx } , , 

■ 2y *1 7 •• P(x>, l x l <a 
^ Y 20 y 


subject to 

/ a ft ( t) dt = 0 
-a T 


( 2 . 1 ) 
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where 


9v, (t , 0) 

<Mt) = , 1 1 1 <a , y = 0 (2.2) 

and v^(x,y) is the displacement component in y direction corre- 
sponding to the layer. Here, Ex, Ey are the Young's Moduli; 

are Poisson's ratios and G v „ is the shear modulus. The 
xy r yx xy 

bielastic constants y^ are given in Appendix A and [1] and the 
kernel k(x,t) is given in Appendix B. The input function p(x) 
is the crack surface pressure and corresponds to the perturbation 
problem. (2.1) is given for the generalized plane stress. The 
plane strain case is obtained by simply interchanging certain 
constants (see [1] ) . 

However, the Fredholm kernel k(x,t) becomes unbounded for 
a = h which represents a broken laminate. It can easily be 
shown that in this case, the part of k(x,t) which contributes 
to further singularities can be expressed as (see Appendix B) 

k (x, t) = /(k*(x,n)e‘ (h “ t)n/5 5' / l a, ll 

s 77 J 20 o 1 

+ k*(x,n)e" (h “ t)n/5 5 /l “3l]a^ (2.3) 

2 

where k?(x,n) and k*(x,rj) are the asymptotic expressions of 

J- ti' 

ki(x,n) and k 2 (x,n) as n + 00 • 

Hence, equations (2.1) can now be written as 


- /?? Mil^t + f a k , f t)<Mt)dt + /jj 
tt -h t-x -n f -h 


f 1 .. V V ) 

- — *y p( X ) 


2y 20 E 1' 


k f (x, t) <J> (t) dt 
|x|<h 


*(t)at = o 
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(2.4) 


where 


k f (x,t) = k(x,t) - k (x, t) = 7 {[k, (x,n) - k*(x,n)] 

f s ttY 20 o 1 1 

e -(h-t)n/55/|»il + [k2(x#n) - k*{x,n)]e'' (h - t)n/ ^5 / l w 3l 

} dn " «> 

It is quite easy to show that k^{x,t) is a bounded kernel and 
thus can be evaluated numerically from (2.5). 

On the other hand, k_(x,t) can be obtained in closed form. 

s 

From the Appendix B, the asymptotic expressions of k 1 (x,n) and 
k 2 (x,n) can be given as follows: 


- 

Iuj, jnh 

- 

| oj 3 1 nh 

k*(x,n) = Y 51 e 

cosh (to nx) 
1 

+ Y 52 e 

cosh (to^nx) 


1 wjnh 

— 

t 

| oj_ | nh 

k* (x,n) = Y 53 e 

cosh (01 Tix) 

+ ^54 e 

cosh(w 2 nx) (2.6) 


Hence, after intermediate manipulations we arrive at 

k s (x ' t} = m c j [ t- (ajh+bjX) + t- (a j h-b j x) 1 


j ~ 1 , . • 4 


(2.7) 


where bielastic constants a j' b j an ^ c j are given in Appendix A. 
In dimensionless variables (2.4) becomes 


1 f} {_!_ + 1 


3 

E c. [ 


•]} 


tt -1 t-x 2y 2q . j = 1 j l T-(a_.+bjX) T ~ (a j“ b j x) 

(j) 0 (T)dT + h/J k f (hx/ht) 4» 0 (t) dx = g(x) r 1x1*1 


/£ <l> Q (T)dT = 0 


( 2 . 8 ) 
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where 


t = t/h , x = x/h 

(1 “ v xv v vx > 

4> d (t) « <Mt) , g(x) = P< h X> ' 

° ZT 2(Ty 


(2.9) 


The solution of (2.8) can be expressed as follows [51; 

»o |T > ° F( 2 > y “ (t/h) . Q [t) , 1 T | <1 (2.10) 

(1-tV (t-h) Y (t+h) Y 

where the unknown function F(t) is Holder continuous in the in- 
terval -1 <t< 1 and o<Re(y)<l. 

The strength of singularity y can be determined by investi- 
gating the dominant part of (2.8) or (2.4) near the end j. ints 
(see [5]). Defining the sectionally holomorphic function 


« (z) . 1 r h tiyat 

* lz ' IT ih t-z 

following results can be obtained: 


(2.11) 


F(-l) _ F ( 1) 


<Mz.) = (~) cotry [■ ■ 1 - ---- - I + *,(z n ) , -h < z =x < h 

1 2 (h+x) Y (h-x) 7 11 X 


h 


4 >(z 2j ) - - (357) 


y- 


F ( 1) 


j sin7ry(h+x) 


h J 


" ( 2E"> 


F (1) 


j sinry(h-x) 


~+$ 2 (z 2j ) , h < z 2j =a.h+b.x < 


(1 + 2b j ) h 


?+ t 3 (z 3j ) , h < z 3 j=a jh-b jX < 


(1 + 2b j ) h (2.12) 

9 — 1/2,3 

where <f> i , i = 1,2,3 are bounded functions everywhere except the 


end points where they behave as follows: 


(z, ) 

*1< Z 1> < I z^hl^ 


D. 


, « 2 (z 2j ) < 


z 2 j“ h ! a 2 ' 


5 




(2.13) 


I z 3j" h l °3 

where ctj and Dj are real constants and ctj<Re(Y). Also, $*(z^) 
satisfies the Holder condition near and at the end points. 

Substitution from (2.12) into the dominant part of either 
(2.4) or (2.8) and considering that F(t) is an odd function, 
the following characteristic equation is obtained to determine y: 

3 c * 

costtY + 7 ~=~ E — =0 (2.14) 

^ Y 20 j=l b . 

D 

which can be shown to have at least one real root. 


After determining y, (2.8) can be solved using the numeri- 
cal method given in [43, i.e. 


i! =1 F( VVi7xI + * hk S < h X i ' hT j> 


■ 9( Xi> 


+ TThk f (hXi/hTj)] 
X ~ 1 , . . ,N-1 


N 

E W.F(t.) = 0 
j=l 3 3 


(2.15) 


where 


p N ( -Y,-Y) ( t _.) » o j = 1, . . ,N 

P N-1 (1 ~ Y/1_Y) (X i 5 = 0 1 = 1 , . • , N-l (2.16) 

and the Wj are the corresponding weights. From (2.15) N unknowns 
p (Tj), j = 1,..,N can be determined. 

3. THE STRESS INTENSITY FACTOR 


The stress intensity factor at the crack tip will be de- 
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fined as follows [ T X y 2 ^ x '^ = : 

K * lim [2(x-h)] Y a (x,0) (3.1) 

x->-h y 2 

x>h 

where a 9 (x,y) corresponds to the matrix. a (x,0) can be ob- 
y ^ Y 2 

tained in terms of <J> (t) starting from the expressions given in 
[1] . Hence, after intermediate steps we obtain, 


^19 (1 - V xy V ?x ) 


dn 


E* 

y 


a y2 (x,0) = <Mt)dt f ^ 


- (h-tjn/^s/lu, | -jws|(x-h)n y?e -|w?|(} 

{e tVl 6 f 13 (n)e + _^f 17 (n)e 


4* e 


-(h-t)n/ 05 /|M 3 | ^ -|w*|(x-h)n y J 5 


[y* f 14 (n)e 


16 


+ Tir f is (,i> 


- | w* | (x-h)n (3.2) 

el ]} 


where 


f 

f 

f 

f 

f 

f 

f 

f 


11 

12 

13 

14 

15 

16 

17 

18 


(n) 

(n) 

<n) 

(n) 

(n) 

(n) 

(n) 

(n) 


^7 

- — + Y 9 o tanh (w. nh) 

'21 * 1 
3o 

Y — + y 22 tan1n ( u 3 r lh) 

2 X 

y 53 f(n) + f ? (n)f 11 (n) + f 8 (n)f 12 (n) 

y 54 f(n) + f 9 (n)f 11 (n) + f io (n) f i2 {r|) 

8 7 + 8^ sign (wS) ' tanh (w^nh) 

6 3 + sign(w|) tanh(w 3 nh) 

Ygif (n) - f 7 (n)f 15 (n) - f g (n) f 16 (n) 
Y 52 f (n) - f 9 (n)f 15 (n) - f io {r| ) f i6 


(3.3) 


(3.4) 


-h) n 

) 
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However# the right hand side of (3.2) becomes unbounded 
near the tip of the crack. Hence, the dominant part which con- 
tributes to the singular behavior can be extracted in a manner 


similar to the procedure used to obtain (2.4). Thus, a (x,0) 

y ** 

can be expressed as: 


(1 ~ W , 

E* 

y 


(x 0) _ f Zl r h 0. ( t) dt . o . 0) 

Y2 (X,0) j =1 tt ih y2 (X ' 0) 


t- (ajh+bjx) 


(3.5) 


where a°^(x,0) is a bounded function. 

To determine the stress intensity factor K, the behavior 
of the Cauchy integral in (3.5) near the end points has to be 
investigated. Following the method given by Muskhelishvili [5] , 
the sectionally holomorphic function defined in (2.11) can be 
expressed near x “ h as follows (h<x<°°) : 


<!> (z . ) 
43 


- (^) Y 


F ( 1) 


Y + V z 4-i } ' 
2b* sirnry (x-h) T - 1 


h<z, .=a?h+b$x < 
43 3 3 


(3.6) 


Here, 4> (z. •) is bounded everywhere except the end points, and 
o 4 3 

near x = h 


♦o (z 4j>l < 1 773 7 

l z 4;f I O 


, a Q <Re(Y) 


where C and a Q are real constants. Hence from (3.1), (3.5) and 


(3.6) we obtain 
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K 


h Y F(l)E* 

* * 

( 1 " v V w v . ) simry 
x y yx 


4 

Z 

j=l 



generalized (3.7) 
plane stress 


For plane strain 
(see [1] ) . 


* * 
1— v v 
xy yx 
* 

E.. 


= E 


x 


•k 


A 


* 

should be replaced by I/A 22 


4. NUMERICAL RESULTS AND CONCLUSION 


The numerical results will be given for p(x) = p Q = constant 
which corresponds to uniform crack surface pressure. H/p^^ 
values as well as y will be given for different material combi- 
nations. 

The following materials will be selected: 

(A) Boron-Epoxy : 

E^ = 3.5 x lO^psi , E y = 3.24 x 10 7 psi 

v =0.23 , G = i. 23 x 10 6 psi 

y A 

and for the plane strain 

E z = 3.5 x 10 6 psi , v zx = v zy = 0.23 

(B) Boron-Epoxy : 

E x = 2.72 x 10 7 psi , E y = 5.5 x 10 6 psi 

v yx = 0.22 , G xy .= 7.0 x 10 5 psi 

and for the plane strain 

E^ = 2.72 x 10 7 psi = V = 0.22 

z zx zy 

(C) Boron-Epoxy : 

E x = 5.5 x 10 6 psi , E y = 2.72 x 10 7 psi 
V yx = 0.1084 , G xy = 7.0 x 10 5 psi 

and for the plane strain 

E,, " 5.5 x lO^psi , y - v = 0.1084 

Z ZX c* jf 

(D) Glass-Fiber (20% volume fraction) : 

E x = 6.6 x 10 5 psi , E y = 2.52 x 10 6 psi 
V yx = 0.32 , G xy = 2.9 x 10 5 psi 

and for the plane strain 

E z = 6.6 x 10 5 psi , u zx = v zy = 0.32 
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Both y and K/p Q h Y are given in Table 1 for various material 
combinations* As seen from the equation (3.7), the stress in- 
tensity factor is determined in terms of F(l) . Since the un- 
knowns obtained from (2.15) are F(tj)(j = 1,..,N), F(l) is found 
by an extrapolation. Hence, three different values for K/p Q h Y 
are given in Table 1 to illustrate the effect of extrapolation. 
Columns 1, 2 and 3 are obtained by extrapolating F(t_.) (j = 1,..,N), 
(j = 1,2,3) and (j = 2,3,4) respectively. From these results, it 
appears that F^) (as well as F(t n )) has a significant effect on 
the stress intensity factor. It is also observed that the y 
and the K/p Q h Y values are quite sensitive to the material properties. 
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APPENDIX a 


For Material Type I s 

Dimensionless bielastic constants y . (i = 46,70); 
(see reference [ 1 ] for the constant^ y^ - y^ 5 ) . 

y 46 ~ [y 27 + y 26 Bign(« x )HY 25 + Y 24 sign(o> 3 )3 

“ (Y 29 + Y 28 sign(a) 3 ) ] [y 23 + y 24 sign ( uj^^) 3 

y 47 = Y 38 + Y 39 P i g n ^ 3 ) 
y 48 = Y 40 + Y 41 sign < w l> 

Y 49 = Y 42 + Y 4 3 ai 9 n < w 3 ) 

Y 50 = y 44 + y 4 5 si 9 n ( w i) 

Y 51 = Y 7 Y 47 /{Y 19' Y 46 ) 

Y 52 = Y 8 Y 48 //{Y 19 Y 46 ) 

Y 53 = Y 7 Y 49 /(Y 19 Y 46 ) 

Y 54 = Y 8 Y 50 //(Y 19 Y 46 ) 

Y 55 = B 9 sign (w^) - 3 * sign(w*) 

Y 11 

Y 56 = y ~ ^8 si< ?n(w$) ~ 0 io sign(w 3 )3 

Y 57 = Y 22 ” B 9 sign([ 0 l ) 

Y 58 = Y^J Y 22 + 6 10 S 1 5 n ^ 3 ' J 

B 7 

Y 59 Y 21 + Y 22 S 1 9 n ^i ) 
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+ y 22 sign(u) 3 ) 


3 S 

Y60 = v 2 i 
Y 61 = ^7 + e g sign (u*) sign{to 1 ) 

Y 62 = ^8 + si 9n(w|) sign (w^) 

Y 63 " Y 57 Y 46 + Y 47 Y 59 + Y 48 Y 60 
Y 64 s Y 58 Y 4 6 + y 49 y 59 + Y 50 Y 60 
Y 65 " Y 55 Y 4 6 " Y 61 Y 47 " Y 62 Y 48 


Y 66 “ y 56 y 46 “ Y 61 Y 49 " Y 62 Y 50 
V 67 = Y 16 Y 63 /{Y 46 Y 19 ) 

Y 68 ” Y 65* Y 15 / ^ Y 46 Y 19 Y 21^ 


y 69 ” y 16 y 64 ^ ^ Y 46 Y 19 ^ 


Y 70 Y 66 Y 15/ (Y 46 Y 19 Y 21 ) 


Bielasttc constants c. , b . , and a. (i = 1,3): 


0) 


C 1 = - 


C 2 = " 


C 3 = ~ 


lL 


w, 2 




Y 51 


H - 1 

, b, = — 


^5 


l w l 1^52+1^,17 


1| ’ 52- r3| y 53 , r i| [ (1) 3 | 


' b 2 ° 




CD 


_ 3 ± y b - C "3 2 

/I? Y54 ' 3 “ ^ 


a . = 1 + d . 
i l 
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I— * 


Bielastic constants c*, b*, a* (i ~ 1,4) 

i 3, 



c 


* 

2 


c 


* 

3 


/37 
I “3 | 



'68 

Y 69 


70 




AP?ENDIX B 


Fredholm kernel (Bielastic constants Y^ are defined in [ 1 J ) 

, 00 ~(h“t)n/B7/| u, | 

k(x,t) = — — / [k,(x,n)e 
vy 20 o 

- (h-t)n/0T/|w 3 | 

+ k 2 (x,n) e ] dn 

k i (x,n) s 2 ' y" I T n T fY 7 f 5 (x*n)f 7 (n) + Y 8 f 6 (x#n> f 8 (n) J 
19 

k 2 (x,n) = 2 y £ (n V ^ Y 7 f 8 n ) f 9 ( n ) + Y 8 ^g (x,ii)£^q (n) 1 

1 19 

f(n) = f 3 (n)f 1 (n) - f 2 (n)f 4 (n) 

f i(n) = Y 2 5 + Y 2 4 tanh (“ghh) 
f 2 (n ) “ Y 29 + Y 28 tanh(w 3 nh) 
f 3 (n) = Y 27 + Y 26 tanh(ai 1 nh) 

f 4 (n ) - Y 23 + Y 24 tanh (w^nh) 
f^fH/X) = cosh (tutrix) /cosh (w^nh) 
fg(ri/X) = cosh (co 3 nx) /cosh (w 3 nh) 
f 7 (n ) = Y 38 + Y 39 ta'nh (w 3 nh) 
f 8 (n) = Y 40 + Y 41 tanh (Wj^nh) 
f 9 (n) = Y 42 + Y 43 tanh(w 3 nh) 
f 10 (H) = Y 44 + Y 45 tanh(w 1 nh) 
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